s— wave superconductivity in the presence of local Coulomb correlations 
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We derive the superconductiong mean-field equations for an attractive interaction, V, in the 
s-wave channel when local Coulomb interactions are taken into account for any value of U. Our 
results show that the Coulomb repulsion is detrimental to the critical temperature, T c , and the 
order parameter, A(T), for values of U > \V\. Furthermore, our results depend on band filling in 
a sensible way. In the presence of local correlations, 2A(0)/T C differs from the BCS ratio, since 
Coulomb interactions affect much more the superconducting critical temperature, T c , than the 
superconducting order parameter, A(T). We conclude that the presence of Coulomb interactions 
play an additional role in the analysis of experimental data, specially in narrow band systems. 
Pacs numbers: 74.20.Fg, 74.60.-w, 74.72.-h 



The discovery of high-T c superconductors (HTSC) [|J has given a huge impetus to the theory of correlation effects 
H due to the fact that there is the belief |j] that the normal properties of these materials could be explained in 
the framework of the Hubbard model since electron correlations are strong, i.e., the on-site electron-electron 

repulsions U are much larger than the energies associated to the hybridization of atomic orbitals belonging to different 
atoms Two dimensional (2 — D) Mott-Hubbard insulators exhibit unconventional electronic, optical, and magnetic 
behavior when doped with mobile charge carriers. However, subsequent analytical and numerical work lead to the 
conclusion that the possibility of superconductivity order out of purely repulsive interactions is not conclusive. So, 
the mechanism of the HTSC is still elusive after almost thirteen years of intense research and extensions to the pure 
Hubbard model have been used 0. See also Ref. ||. Due to this lack of consensus an additional phenomenological 
interaction has to be included to go into the superconducting phase. Because of that, following Ref. \7\, we postulate a 
Hamiltonian which is composed of two terms, the Hubbard contribution, which we call Hjj and the Cooper part which 
we represent for Hy. To say the truth, this Hamiltonian should be considered the minimum model to cover a great part 
of the rich phase diagram of the cuprate superconductors which include antiferromagnetism (AF) , superconductivity 
(SC), insulating state (I), non-Fermi liquid behavior (NFL), lattice distortion, etc. ||. 

Due to the fact that there has been a large amount of studies in the normal state of the HTSC, we consider that the 
normal state Green function is known and we tackle the superconducting state by means of a mean-field treatment. 
For the normal state Green function, for all values of U , we choose an academic Hubbard-III like approximation JTcj ] 
which gives a metal for small values of U (no gap in the density of states) and becomes an insulator for large values 
of U (there is a gap for U > W, where W is the band width). Our studies differ from the ones recently published 
in the literature where mean field analysis has been performed for both U and V Jn],[l^] . We then find the diagonal 
and off-diagonal superconducting one-particle Green functions, G(k,iuj n ) and F(k,iuj n ), respectively. Using these 
two Green functions, we derive the self-consistent equations for the density and the gap equation, respectively. These 
equations are self-consistent and we have solved them numerically. 

In the superconducting phase, we will follow the BCS treatment supposing that the presence of correlations leave 
this formalism untouched S. Our dynamical equations become 

G^ 1 (k, uu n )G(k, nu n ) + A(T)Ft (k, iLu n ) = 1 , 

Gv l (%,iu n )F*(%,iu n )- A{T)G(k,iLu n ) = , (f) 
with Gu(k,iuj n ) being the normal state correlated one-particle Green function |UJ.M 



Gu(k,iu; n ) = - 1 9 - + - 9 tj , (2) 

vjj n + fi — e(k) iui n + fi — e(k) - U 

e(k) being the free band structure, p is the particle density/spin and u> is the fermionic (odd) Matsubara frequency. 
According to our interpretation, the Hamiltonian is split in two pieces 
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H = Hr 



(3) 



where Hu — H a is our unperturbed Hamiltonian for which we know the solution in the normal state. This is the 
reason of having chosen the normal Green function as in Eq. (|^). Of course, we do not know the exact solution of the 
Hubbard model even in d = oo due to the local character of the Coulomb interaction However, if we have a normal 
Green function which interpolates between the weak and the strong coupling regime, then we can say that we have a 
pretty good solution to the the Hubbard part. According to this view what we are studying is the effect of Coulomb 
interactions, U, on superconductivity. We mention that there are other types of approximations for Gjj(k,iuj n ) as it 
has been done in Ref. Oh i.e., Hartree-Fock (HF), second order perturbation theory (SOPT) and the alloy analogy 
approximation (AAA) pi]. The way to view Eq. (||) in terms of Feynman diagrams is the following: our internal 
lines are due to correlated Green functions, due to Gu(k,iu n ), and the vertex or interaction is due to the Cooper 
interaction, V. 

Solving Eqs. (fy), we get 
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Gy 1 (fc,^„)G^ 1 (fc, -iu n ) + \A(T)\ 
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G^(k,iu, n )G^(k, -iw n ) + |A(T)| 
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Before showing our working equations, it is instructive to elaborate a little further our superconducting Green 
functions, expressing them in terms of their poles and their residues. The result of this is the following 
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The notation is explained in Eqs. (j^). The reason of having four poles is because the denominator of G(fc, iuj n ) and 
F'(k, iu> n ) is a polinomial of order four. 

Substituting Gu{k,iui n ) from Eq. (0) we find the mean-field equations for A(T) and p, respectively, as 
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P = TTFi / MaO/(Yi(aO) + a 2 (x)f(-Yi(x)) + a 3 (x)f(Y 2 (x)) + a A {x) f (-Y 2 {x))) dx 
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where we have chosen a flat free density of states, i.e., i.e., Nl(c) = 1/2D for — D < e < +D and zero otherwise. In 
the end we have chosen 2D = 1. The notation in Eq. (||) is as follows 
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From our previous expressions we can verify that J_5| 

oti + a 2 + a 3 + o;4 = 1 ; Yi(ai - a 2 ) + Y 2 (a 3 - a A ) = e(k) + pU 



(8) 



Eqs. (JsJ) are nothing that the first two sum rules or moments of the diagonal spectral function, A(k, lu) = 
T/7rlim ( 5^Q+ Im G(k,oj + iS) . We see that the presence of Coulomb correlations yields a\ ^ a 2 and a 3 ^ a 4 . 
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Analyzing our gap equation we see that the BCS energy symmetry is kept, i.e., we have the following one-particle 
energy excitation solutions, ±Yi, and ±Yz (the last line of Eq. (0)). By the same token, the off-diagonal spectral 

function, B(k, u>) = — 1/tt lim,5_ >0 + Im F(k, uj + iS) complies with the first two off-diagonal sum rules. These consid- 
erations are a check that we are working in a mean field scheme for the superconducting one-particle Green functions, 
as it should be. To go beyond the mean-field solution, we need to include pairing fluctuations as it has been done by 
Micnas et al of Ref. Also, pairing fluctuations have been studied by Schmid |L6| and others Jl7| , ^8 | 

In Fig. 1 we present T c vs U for various values of fi, i.e., fi = 0.25; 0.50; 0.75 and 1.00, for several values of V, i.e., 
V = —0.50; —1.00; —1.50 and —2.00. We inmediately conclude that T c goes to zero for high values of U. Also, we 
observe that, for chosen values of V, there is a maximum value of U ~ 1.4 beyond which there is not superconductivity. 
We gain these results by taking A(T) = in Eqs. ([!]). Fig. 1 is our guiding line for solving Eqs.(|l|) below T c , i.e., we 
will use the parameters of Fig. 1 and solve A(T) ^ for T < T c . 

For example, in Fig. 2 we plot A(T)/U vs T/T c for various values of U and V for p = 0.50. For all our curves, 
A(T)/U goes down when U increases. This can be explained because A(T) is a parameter which is conceptually 
defined thru \V\, i.e., it depends little on U. When we normalize A(T) with U, we are decreasing this ratio drastically. 
For high values of U, A(T)/U goes to zero. We see that for U = 1.0 the value of the normalized gap, i.e., A(T)/U, 
is almost zero. However, if the value of \V\ is comparable to U then we can have a sizeable value of the normalized 
ratio. So, there are two competing parameters in the theory, U and V, which are controling the value of the order 
parameter. 

In Fig. 3 we present A(T)/A max vs T/T c , where A max = A(0). T c is the value given in Fig. 1. We see that 
the outmost curve is the one with U = (pure BCS case) and the inmost curve is the one with the highest value 
of U . Thus, for high values of U, A(T)/A max decreases in the intermediate region, i.e., for < T < T c . So, 
superconductivity is basically diminished for high values of U, a result which had been reached in previous works 
|^l|,Q. In particular, for p — 0.50, and U = —1.00 there is a great deviation with respect to the pure BCS case 
(U = 0). This fact should be detected experimentally. Fig. 3 is particularly iluminating because depending on band 
filling (or p) and certain values of U we can deviate from the pure BCS result. Of course, there is an additional 
difference with respect to the U = case: we have two different excitation energies, namely Y± and Yi (Eq. (^)) 

In addition, we have evaluated 2A(0)/T C . In pure BCS, it is approximately equal to 3.5. For p = 0.25, V = —0.50, 
U = 0.25 (p w 0.6, we get that this ratio is w 4.37. Also, for p = 0.50, V = -0.50, U = 0.75 (p w 0.57) this ratio is 
~ 6.585. As we have mentioned before, A(T) is a quantity which is defined mainly thru the pairing interaction, V. 
So, it does not change too much with U. On the contrary, U has a bigger effect on T c reducing it. This argument 
may explain the reasons why the 2A(0)/T C can be very different from the pure BCS ratio. In other words, the 
presence of Coulomb interactions plays a fundamental role in changing the universal ratio. We also see that the 
band filling, in presence of Coulomb interactions, plays a sensible role too. These results may have some relationship 
with one of the characteristic features of superconductivity in the cuprate superconductors: the strong dependence 
of T c and A(T) on doping concentration, x, away from half filling (n = 2p = 1 — x) JljJ. But, above all, what 
we are led to is to see that the ratio 2A(0)/T C is not longer universal and changes with doping, for example. This 
has very important experimental consequences since we do not need to consider non-BCS* |18| theories to explain 
these big ratios. However, we can say that in our approach the Coulomb interaction plays an equivalent role of the 
Non Fermi Liquid parameter a of Refs. |l8[ , namely, they decrease the value of T c . 

In short, we have solved the gap and density equations inside the superconducting phase for s-wave symmetry order 
parameter taking into account local Coulomb repulsion of any strength. Speaking a little bit about numerics, we say 
that we have paid due care to the cases Yi = Y2 and Y2 — 0. Also, for small values of A(T) convergence depends on 
the initial conditions. The presence of electron correlations is detrimental to superconductivity When we say that 
U conspire against superconductivity, one is saying that it diminishes both T c and A(T). For the first part of this 
statement, we show Fig. 1, which we use as our reference. The second part of the statement is shown in Figs. 2 and 
3. At this point we call the attention of the reader to the point that our calculations were done at fixed \i. When U 
is taken into account, n (or band filling) starts to play a decisive role. In particular, we have found 2A(0)/T C > 3.5 
which has also been found in experimental measurements. We could tried another type of normal state one-particle 
Green functions, instead of the one of Eq. (||). We leave this task for the future. We have said that our choice for 
the normal one-particle Green function, G]\[(k, iu> n ), is an academic one because the weights of the spectral function, 
A(k,u), are 1 — p and p, respectively. In other words, we do not have fc-dependence which is a vital ingredient to 
explain angle-resolved photoemission (ARPES) data |2(| of HTSC. However, our aim was to study the effect of U 
on the superconductivity parameters T c and A(T). Then, other better approximations can be used as it has been 
done in Ref. |ll|. Furthermore, one could use the two-pole Ansatz of Nolting |ll| which may give a normal state 
metal-insulator transition if the band narrowing band factor, B s {k), is properly treated. However, the basis physics 
will remain the same. Another point which has not been considered here is the effect of the lattice structure El]] 
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(nearest-neighbor (n.n.), second-nearest-neighbor (n.n.n.), etc. A recent use of the two-pole ansatz for high values 
of U to interpret exact diagonalization data of two dimensional clusters is given in Ref. |g2| . The effect of U on the 
isotope coefficient is a problem worthwhile to explore p3] as well as the jump on the specific heat at T = T c [^H . 
This paper has dealed with the s-wave order symmetry only. However, complex symmetry [^5| order parameter |25f| 
do not represent numerical difficulty. See also, Ref. p6| . 
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Figure Captions 

Figure 1. T c vs U, for different values of V, i.e., V = -0.50; -1.00; -1.50 and -2.00. (a) fi = 0.25; (b) fi = 0.50; (c) 
fi = 0.75; and (d) /z = 1.00. 

Figure 2. A{T)/U vs T/T c for several values of U and four values of V, i.e., (a) V = -0.50; (b) V = -1.00; (c) 

V = — 1.50; (d) V = — 2.00. Here we have fixed the chemical potential to fi = 0.50 

Figure 3. A(T)/A max vs T/T c for several values of U and four values of V, i.e., (a) V = -0.50; (b) V = -1.00; (c) 

V = —1.50; (d) V — —2.00. We have again fixed the chemical potential to fj, = 0.25 
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